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Resumen 


La ganadería en el piedemonte llanero colombiano está ampliamente distribuida y se desconoce su impacto sobre 
los murciélagos insectívoros aéreos. Monitoreamos métodos acüsticos pasivos a estos murciélagos en un sistema 
silvopastoril (SSP) plantado con Mimosa trianae, un sistema ganadero convencional dominado por pasto (SC) y 
un bosque (B) en el piedemonte llanero colombiano. La composición de especies del B se diferenció de las otras 
coberturas, mientras que no hubo diferencias entre SC y SSP. Hubo mayor riqueza en el SC que en las otras co- 
berturas, mientras que no hubo diferencias entre el B y el SSP. Myotis nigricans, Peropteryx macrotis, Saccopteryx 
leptura, Molossus cf. rufus, Eptesicus I y Cormura brevirostris tuvieron la mayor actividad en el paisaje, y su uso 
del espacio dependió de la cobertura. La mayoría de las especies usaron el SC que las otras coberturas. En con- 
clusión, la heterogeneidad espacial asociada al paisaje ganadero estudiado tiene efectos tanto en la estructura 
del ensamble como en las respuestas de las especies que lo conforman. 


Palabras clave. Agroecosistemas. Bioacústica. Chiroptera. Orinoquia. Sistema silvopastoril. 


Abstract 


Cattle ranching is widely spread in the foothills of the Eastern Cordillera of Colombia, and its impact on aeri- 
al insectivorous bats is unknown. We used passive-acoustic monitoring of these bats in a silvopastoral system 
(SSP) planted with Mimosa trianae, a conventional livestock system dominated by pasture (SC) and a forest (B) in 
the Eastern Cordillera foothill of Colombia. Species composition of B differed from SC and SSP, whereas there 
were no differences between SC and SSP. Species richness was higher in the SC than in the other land covers, 
whereas there were no differences between B and SSP. Myotis nigricans, Peropteryx macrotis, Saccopteryx leptura, 
Molossus cf. rufus, Eptesicus I and Cormura brevirostris had the highest activity in the landscape, and their habi- 
tat use depended on land cover. Most bat species used SC more than the other land covers. We found that the 
spatial heterogeneity associated with the studied cattle ranching landscape affects both the structure of the bat 
ensemble and the responses of its species. 
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Introducción 


La ganadería es una de las principales actividades cau- 
santes de la transformación del paisaje en el planeta 
(Ramankutty ef al., 2006), y se estima que aproximada- 
mente un tercio de la superficie terrestre del planeta está 
destinada para la ganadería y la agricultura (Ellis et al., 
2010). En Colombia la situación es similar, y se ha estima- 
do que un 30 % del territorio está dedicado a la ganade- 
ría, aproximadamente 34 millones de hectáreas, y más del 
77 % de esa área sufre algún grado de erosión (IDEAM 
& UDCA, 2015). Así, la expansión ganadera genera una 
transformación de las coberturas naturales, incrementa 
la presión sobre los recursos naturales y tiene un efecto 
negativo sobre la biodiversidad (Alkemade et al., 2013; 
Dettenmaier ef al., 2017). La ganadería está entre las pri- 
meras causas de pérdida de biodiversidad en Colombia, 
y en el piedemonte llanero amenaza los ecosistemas de 
bosque (Romero et al., 2008; Etter et al., 2017). Por todo 
lo anterior, se ha resaltado la necesidad de estrategias 
para conservar la alta diversidad del país, aun en esos 
ambientes altamente perturbados (Moreno et al., 2018). 


Para contrarrestar los efectos negativos de la ganade- 
ría convencional, extensiva y basada principalmente 
en el uso de pastos exóticos, se ha sugerido el uso de 
los sistemas silvopastoriles, SSP (Calle et al., 2012), que 
integran el uso de árboles o arbustos con pastos para el 
ganado (Murgueitio ef al., 2013). En comparación con 
los sistemas convencionales de ganadería, la cobertura 
heterogénea de los SSP puede facilitar la recuperación y 
mejoramiento de suelos, al mejorar el flujo de nutrientes 
y los ciclos locales del agua, (Alonso, 2011; Naranjo et al., 
2012). En términos de los efectos sobre la biodiversidad, 
los arbustos y cercas vivas de los SSP pueden funcionar 
como corredores que favorecen la presencia de especies 
de murciélagos, aves, escarabajos, mariposas y fauna del 
subsuelo (Harvey et al., 2005; Mahecha, 2002). 


En Colombia, los estudios sobre murciélagos en áreas 
ganaderas se han realizado en el Caribe, y se ha encon- 
trado una mayor diversidad de Phyllostomidae en los 
SSP que en los sistemas convencionales (Ballesteros, 
2015). También, una especie de murciélago frugívo- 
ro ha mostrado mejor condición corporal en los SSP 
que en sistemas convencionales, y algunas especies 
pueden acumular metales pesados producto de la adi- 
ción de agroquímicos a los SSP (Chacón-Pacheco & 
Ballesteros-Correa, 2019; Racero-Casarrubia et al., 2017). 
Así, es limitada la información sobre murciélagos in- 
sectívoros aéreos, 1. e., que se alimentan de insectos en 
el aire (Wilson, 1979). 


Murciélagos insectívoros aéreos en un paisaje ganadero 


Algunos murciélagos insectívoros aéreos toleran am- 
bientes altamente modificados y allí ayudan a conser- 
var funciones ecosistémicas e incluso pueden prestar 
servicios ecosistémicos al controlar insectos que afec- 
tan negativamente los sistemas de producción (Boyles 
et al., 2011; Cleveland et al., 2006). Entre las herramientas 
para permitir que los murciélagos insectívoros persistan 
en ambientes rurales se encuentra el manejo de áreas 
arboladas (Guldin et al., 2007). Esto se debe a que las 
zonas arboladas naturales en paisajes altamente mo- 
dificados por los humanos pueden afectar el uso del 
espacio de los murciélagos insectívoros aéreos, dado 
que pueden servir como refugio, fuentes de alimento 
o como corredores de tránsito (Barclay & Kurta, 2007; 
Grindal & Brigham, 1999). Además, en ambientes alta- 
mente intervenidos, los murciélagos hacen uso diferen- 
cial del espacio, segün las coberturas presentes, ya que 
estas pueden afectar la disponibilidad de presas y las 
oportunidades de forrajeo (Avila-Flores & Fenton, 2005; 
Jung & Kalko, 2010). Así en diferentes localidades del 
neotrópico, se ha encontrado que cambios en el paisaje 
producto de la actividad agropecuaria pueden afectar el 
uso del espacio por los murciélagos insectívoros aéreos 
(Estrada et al., 2006; Williams-Guillén & Perfecto, 2011). 


Las comunidades de murciélagos en la base de los 
Andes orientales de Colombia se caracterizan por te- 
ner un rico ensamble de especies insectívoras aéreas 
(Ramírez-chaves et al., 2013; Sánchez, 2017); ensamble 
(Ramírez & Gutiérrez-Fonseca, 2016) corresponde a un 
grupo filogenéticamente relacionado, i. e., Chiroptera, 
que explota un recurso similar, i. e., insectos aéreos. En- 
tender cómo los ensambles de murciélagos insectívoros 
aéreos responden a las prácticas ganaderas puede ser- 
vir para adoptar estrategias de manejo que reduzcan el 
impacto de esta actividad económica. Por ello, estudia- 
mos un ensamble de murciélagos insectívoros aéreos 
en el piedemonte llanero colombiano, comparando un 
sistema ganadero convencional (SC), i. e., pastizal, un 
sistema silvopastoril (SSP) dominado por Mimosa tria- 
nae, y un remanente de bosque nativo. Nuestra primera 
hipótesis es que la riqueza y composición de especies de 
murciélagos insectívoros difieren entre las coberturas. 
Esto se debería a que el incremento en la densidad de 
obstáculos asociados a la cantidad de árboles presenta 
para los murciélagos problemas mecánicos por la nece- 
sidad de maniobrar, y problemas de percepción por la 
necesidad de detectar y evitar los obstáculos (Fenton, 
1990). En consecuencia, la densidad de obstáculos en un 
ambiente se relaciona con la variación en la morfología 
de las alas y de las señales de ecolocalización usadas por 
murciélagos insectívoros, ocupando diferentes espacios 
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(Aldridge & Rautenbach, 1987; Thies et al., 1998). Nues- 
tra segunda hipótesis es que los murciélagos insectivo- 
ros aéreos hacen un uso diferencial de las coberturas. 
Esto debido a las razones explicadas anteriormente, y 
porque en México se encontró una mayor actividad 
de murciélagos en fragmentos de bosques que en pas- 
tizales (Estrada et al., 2004), mientras que en Estados 
Unidos se halló una relación positiva entre la actividad 
y la distancia entre árboles (Gerth & Chelsvig, 2008). 


Materiales y métodos 


Área de estudio. El piedemonte llanero colombiano 
es una subregión de la Orinoquia que se ubica en la 
base de la Cordillera Oriental de los Andes (Rangel-Ch. 
et al., 2011) y ha sido intensamente transformada por la 
actividad humana. Realizamos este estudio en la Fin- 
ca La Rosania, ubicada en el municipio de San Luis de 
Cubarral, departamento del Meta (3° 49 25.84" N 73° 
50° 2.87" W, 623 m s.n.m, Figura 1). El municipio tiene 
una precipitación anual de 2000-4000 mm y una tempe- 
ratura anual de 20-27 °C (Callejas Avila, 2017: González 
Rey, 2017). A lo largo del año la precipitación presenta 
un régimen unimodal biestacional, con bajas precipita- 
ciones entre diciembre y febrero, y altas precipitaciones 
entre marzo y noviembre (Datos IDEAM de 10 años, 
estación: 32060030). 


La ganadería y la producción de leche son las principa- 
les actividades productivas de la finca. En el momento 
del estudio, los pastizales eran dominados por el género 
Brachiaria, cubrían aproximadamente 32.5 ha, y tenían 
árboles dispersos con alturas menores a 7 m. El SSP esta- 
ba conformado por árboles de falso yopo Mimosa trianae, 
tenía un área de 1.23 ha, una edad de siembra de 2 a3 
años, y los árboles estaban sembrados cada 5 metros y 
tenían hasta 5 m de altura. Los pastos en el SSP eran 
del género Brachiaria, y allí se removían las ramas secas 
del suelo, para que no interfirieran el desplazamiento 
del ganado. El bosque de piedemonte de la finca tenía 
un área de 106.7 ha, con un estrato arbóreo de 12-25 m 
de altura, un estrato de arbolitos de 5-12 m de altura y 
un estrato arbustivo de 1.5-5 m de altura, según la des- 
cripción de Rangel € Velásquez (1997). 


Muestreo acústico. Realizamos en 2017 un muestreo 
estratificado, con nueve noches en febrero (baja precipi- 
tación), diez noches en junio, cuatro noches en octubre 
y una noche en noviembre (alta precipitación). En total, 
obtuvimos 24 noches efectivas de grabación (288 horas), 
i. e., 8 noches efectivas y 96 horas por cobertura. Usamos 
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un detector ultrasónico SM4BAT FS con un micrófono 
externo SMM-U1 (Wildlife Acoustics® Maynard, MA, 
USA), El detector se instaló en seis sitios de grabación 
por cada cobertura. En cada punto de muestreo ubica- 
mos el detector con el micrófono a 3 m sobre el nivel 
del suelo, en un ángulo de 45° respecto a la superficie 
del suelo y apuntando hacia el interior de la cobertu- 
ra, desde las 18:00 h hasta las 06:00 h del día siguien- 
te. Configuramos el detector para registrar frecuencias 
entre 16 y 120 kHz, con una tasa de muestreo de 256 
kHz, una ganancia de 12 dB y una longitud máxima 
de archivo de grabación de 15 segundos. El detector 
estuvo activo todo el tiempo de instalación, pero solo 
grabó cuando registró señales emitidas en el rango de 
frecuencias programado. Los archivos se guardaron en 
formato WAV sin compresión. 


Dado que no existen librerías de referencia para la Ori- 
noquia colombiana, realizamos identificaciones manua- 
les siguiendo claves acústicas y literatura (Jung et al., 
2014; Arias-A guilar et al., 2018; López-Baucells et al., 
2016; Zamora-Gutiérrez et al., 2016; Zurc et al., 2017) 
y comparando los siguientes parámetros: forma del 
pulso de ecolocalización, alternancia, frecuencia ini- 
cial, final, máxima, mínima, de máxima energía, ampli- 
tud, duración e intervalo entre pulsos. Analizamos las 
grabaciones en el programa BatSound-Sound Analy- 
sis V.4.2., y para las mediciones usamos el programa 
Raven Pro V.1.6 (The Cornell Lab of Ornithology). En 
ambos programas usamos una transformada rápida de 
Fourier (“fast Fourier transform”) de 512 muestras con 
una ventana tipo Hamming. Dado que la configuración 
visual del programa afecta la percepcion de los pulsos 
de ecolocalización, usamos los mismos valores en los 
parámetros para todas las grabaciones analizadas: 1500 
ms por imagen, threshold = 10, amplitud de contraste 
= 3, en BatSound; ajuste de brillo = 63, ajuste de con- 
traste = 67 y overlap = 95% en Raven Pro. Los pulsos de 
ecolocalización no asignados a una especie se registra- 
ron como sonotipos, 1. e., un grupo de sonidos con las 
mismas características y que asumimos correspondían 
a la misma especie (Bader et al., 2015). Los registros con 
baja intensidad o con pulsos muy tenues (< 45 dB) se 
clasificaron como “no identificados”. 


Los murciélagos insectívoros aéreos emiten señales de 
ecolocalización que pueden ser agrupadas en tres fa- 
ses: búsqueda, aproximación y terminal (Schnitzler & 
Kalko, 2001). La emisión de fases de búsqueda repre- 
senta el desplazamiento o la búsqueda de presas, las 
fases de aproximación representan la detección de los 
objetos o presas, y la fase terminal (“feeding buzz”) 
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Figura 1. Ubicación del detector (puntos blancos) para el muestreo acústico de la actividad de murciélagos insectívoros aéreos 


en tres coberturas de un paisaje ganadero del municipio de Cubarral, departamento del Meta, Colombia. 


Figure 1. Location of the detector (white dots) for the acoustic sampling of the activity of aerial insectivorous bats in three land covers of a 


livestock landscape in Cubarral, Meta, Colombia. 


indica intentos de los murciélagos para capturar presas. 
Por ello, evaluamos el uso de las coberturas contando el 
número de fases de búsqueda y las fases terminales de 
cada especie por noche, como un indicador de la activi- 
dad del murciélago en cada cobertura (Sánchez, 2011). 


Análisis estadísticos. Para soportar las identificacio- 
nes asignadas, realizamos análisis de función discrimi- 
nante (DFA) (Parsons & Jones, 2000) por cada familia, 
con base en diez grabaciones de cada identificación y 
medimos tres pulsos por secuencia, de los cuales se ex- 
trajeron los parámetros antes mencionados (Anexo 1). 
Medimos el total de archivos paras las identificaciones 
que no superaron las diez grabaciones y seleccionamos 
grabaciones de las diferentes coberturas. También me- 
dimos los pulsos bajos, medios y altos para las especies 


que presentan alternancia en los pulsos de ecolocali- 
zación. De acuerdo a lo sugerido por López-Baucells 
et al., (2016), agrupamos en Eptesicus I las posibles es- 
pecies de este género. 


Para determinar la variación de la composición de mur- 
ciélagos entre las coberturas e identificar, las especies y 
los sonotipos causantes de la variación, usamos un aná- 
lisis de similaridad ANOSIM, con 9999 permutaciones 
y un análisis SIMPER (Clarke, 1993). Estos análisis se 
realizaron en el programa Past V.3. También construi- 
mos curvas de acumulación de especies con la función 
“specaccum” del paquete “vegan”, para evaluar posi- 
bles diferencias en la riqueza de especies entre cobertu- 
ras, en el programa R versión 3.4 (Oksanen et al., 2019; 
R Core Team, 2019). 
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Evaluamos el uso de coberturas de los murciélagos in- 
sectivoros con un modelo lineal generalizado (GLM) con 
una distribución de Poisson (Sheskin, 2004; Zar, 2010). 
Usamos esta opción, dada la naturaleza de conteo de 
los datos. En el modelo se incluyó el número de fases 
de búsqueda de todas especies registradas como varia- 
ble dependiente, mientras que el tipo de cobertura y el 
mes entraron como variables independientes. También, 
realizamos el mismo GLM usando como variable depen- 
diente el número de fases terminales de todas las especies 
registradas. Para detectar diferencias entre coberturas y 
meses, realizamos contrastes con un ajuste de Bonferro- 
ni (Sheskin, 2004). Debido a la cercanía entre los sitios 
de muestreo, aplicamos la prueba de autocorrelación 
espacial con el índice de Moran, usando una matriz de 
pesos con las distancias entre los cinco vecinos más cer- 
canos. El análisis se realizó con la función “moran.test” 
del parquete “spdep” de R, usando los cinco vecinos más 
cercanos de cada punto de muestreo (Bivand & Wong, 
2013; R Core Team, 2019). Evaluamos, además, el uso del 
espacio para las especies con mayor número de fases de 
búsqueda y terminales, 1. e., especies con más de 150 fases 
de búsqueda en todo el muestreo. Para ello, realizamos 
un GLM con una distribución de Poisson con la misma 
estructura del análisis general del uso de coberturas. 


Resultados 


Análisis del ensamblaje. En total analizamos 10 845 
archivos de audio, correspondientes a 28.8 GB de in- 
formación, y con 6387 archivos con pulsos de ecolo- 
calización de murciélagos. Encontramos 23 sonotipos 
representadas en cuatro familias, Además, siete iden- 
tificaciones a nivel de sonotipo y algunos registros de 
baja intensidad, que no pudieron ser identificados 
(Anexo 2). Para la familia Molossidae, identificamos 
seis sonotipos, dos a nivel de familia, dos a nivel de 
género y dos a nivel de especie. Para Emballonuridae, 
identificamos un sonotipo a nivel de género y cinco 
a nivel de especie. Para la familia Vespertilionidae, 
identificamos seis sonotipos, dos a nivel de familia, 
dos a nivel de género y dos a nivel de especie. Para 
Thyropteridae, identificamos un sonotipo a nivel de 
género. En el SC, registramos 23 especies, de las cuales 
14 fueron exclusivas de esa cobertura, mientras que 
en el B hubo 12 especies, seis exclusivas allí. Para el 
SSP registramos 15 especies, todas compartidas con 
el Bo el SC. En consecuencia, la composición de mur- 
ciélagos de acuerdo con el ANOSIM varió significati- 
vamente entre las coberturas (R = 0.234, p < 0.01). La 
composición del B fue significativamente diferente a 
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la de las otras dos coberturas (B vs. SC: p < 0.001; B vs. 
SSP: p = 0.042), mientras que esto no ocurrió entre el 
SSP y el SC (p = 0.058). El SIMPER indicó que entre el 
B y el SC hubo una disimilitud del 74.9 %, entre el B 
y el SSP, del 72.7 % y entre el SSP y el SC, del 52.7 %. 
Las especies que aportaron en mayor proporción a 
la disimilitud entre las coberturas fueron: Eptesicus I, 
Myotis nigricans, Molossus cf. rufus, Peropteryx macrotis, 
Cormura brevirostris, Saccopteryx leptura (Anexo 3). Por 
otra parte, la riqueza de especies fue mayor en el SC 
que en las otras coberturas, mientras que no hubo di- 
ferencias entre el B y el SSP (Figura 2). 


Análisis general del uso de las coberturas. No hubo 
efecto de la cercanía de los sitios de grabación en la 
frecuencia de las fases de búsqueda y terminales (des- 
viación estándar [d.e] = -0.253; p = 0.600; d.e = -0.338; 
p = 0.368, respectivamente). Por ello no se incluyó este 
factor en los análisis subsecuentes. La frecuencia de fa- 
ses de búsqueda y terminales fue afectada por el tipo de 
cobertura (2.1 = 2; d.e = 999.683; p < 0.001; g.l = 2; d.e = 
73.480; p < 0.001, respectivamente) (Figura 3), siendo 
más usado el SC (x= 395.9 fases de búsqueda; x = 376 
fases terminales), seguido del SSP (x = 232.7 fases de 
búsqueda x= 30.1 fases terminales) y el B (x = 139.2 fases 
búsqueda; x = 15.7 fases terminales). También, la fre- 
cuencia de las fases de búsqueda y terminales dependió 
de los meses (2.1 = 3; d.e = 275.196; p < 0.001; g.l = 3; d.e 
= 231.158; p < 0.001, respectivamente) (Figura 3), sien- 
do mayores los valores en junio (x= 294.7 fases de bús- 
queda; x = 52.2 fases terminales), seguidos de octubre 
(x= 219.0 fases de búsqueda; x = 22.9 fases terminales) 
y por último los de febrero (x = 198.8 fases búsqueda; 
x = 15.0 fases terminales). 


Uso de las coberturas por especie. Todas las espe- 
cies examinadas, dos de la familia Vespertilionidae 
(Eptesicus I, M. nigricans), tres Emballonuridae (P. 
macrotris, S. leptura y C. brevirostris), y un Molossidae 
(M. cf. rufus), presentaron diferencias significativas en 
el uso de las tres coberturas, en términos del número 
de fases de búsqueda y terminales (en todos los MLG's 
p < 0.01). En todos los casos fue mayor el uso del SC, 
seguido por el SSP y por el B (Figura 4), excepto para S. 
leptura. Para este embalonúrido, de acuerdo con las fases 
de búsqueda, el uso fue SSP>SC=B, mientras que, según 
las fases terminales, fue SSP>SC>B. Por otra parte, el 
mes de muestreo también tuvo un efecto significativo 
sobre el número de fases de búsqueda y fases termi- 
nales en las seis especies (en todos los MLG's p < 0.05; 
Figura 4), con excepción del número de fases terminales 
de M. nigricans (2.1 = 3; d.e = 7.332; p > 0.05) (Anexo 4). 
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Figura 2. Curvas de acumulación de especies de murciélagos insectívoros aéreos de tres tipos de coberturas en un paisaje 


ganadero en Cubarral, Meta, Colombia. Las horas de grabación se usaron como unidades de muestreo. La línea continua 


representa la riqueza estimada para cada cobertura y el espacio sombreado de cada línea corresponde al intervalo de con- 


fianza del 95 %. 


Figure 2. Species-accumulation curves of aerial insectivorous bats recorded in three types of land cover in a cattle ranching landscape in Cu- 


barral, Meta, Colombia. Recording hours were used as sampling units. The solid line represents the estimated richness for each cover and 


the shaded area corresponds to the 95 % confidence interval. 
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Figura 3. Uso de coberturas en los meses de muestreo por los murciélagos insectívoros aéreos en un paisaje ganadero en 
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Figure 3. Use of land covers types by aerial insectivorous bats on the sampling months in a cattle ranching landscape in Cubarral, Meta, 


Colombia. Land cover use was measured as number of search phases per night (upper panel). Circles are averages and whiskers are one 


standard deviation. 
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Figura 4. Uso de coberturas por seis especies de murciélagos insectívoros aéreos en un paisaje ganadero en Cubarral, Meta, 


Colombia. El uso se midió según el número de fases de búsqueda (panel superior), y terminales (panel inferior), por noche. 


Los círculos son promedios y los bigotes, una desviación estándar. 


Figure 4. Use of land cover type by six aerial insectivorous bats species in a cattle ranching landscape. Use was measured in number of search 


phases (upper panel) and feeding buzzes (lower panel) per night. Circles are averages and whiskers are one standard deviation. 


Discusión 


Análisis del ensamble. De acuerdo a lo esperado, en- 
contramos diferencias en la composición y riqueza de 
especies entre las coberturas del paisaje ganadero. La 
composición encontrada en el bosque, fue significati- 
vamente diferente a la encontrada en el SSP, pero no 
detectamos diferencias entre esta última cobertura y el 
sistema convencional. Así, en términos de composición 
el SSP es más parecido a un SC que al bosque. Nuestros 
resultados del SIMPER, que evalúan la disimilaridad 
teniendo en cuenta la representación de cada especie 
por cobertura (Mumby, 2001), también indican que las 
especies con mayor número de registros son las que 
aportan más a la diferenciación entre las coberturas. La 
mayoría de esas especies tienen marcadas diferencias 
en el uso de cada cobertura, lo que explica la diferen- 
ciación. Además, la riqueza en el SC, fue seguida por 
la del SSP y la del bosque. Así, los árboles espaciados 
y su bajo porte, ~5 m, en el SSP, permiten que espe- 
cies de espacios abiertos usen esta cobertura. Por otra 
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parte, la densidad de obstáculos y altura de los árboles 
en el bosque actúan como un filtro para algunas espe- 
cies de murciélagos insectívoros aéreos, es decir que 
hay un efecto de la estructura de la vegetación sobre 
el ensamblaje de murciélagos insectívoros (Jung et al., 
2012; Marques et al., 2016). En otras palabras, estos re- 
sultados concuerdan con el hecho que hay murciélagos 
insectívoros aéreos que prefieren el interior del bos- 
que, y otros que prefieren los espacios abiertos (Jung 
et al., 2012; Schnitzler & Kalko, 2001), pero a nivel de 
composición no hay una clara diferencia en la cobertu- 
ra con una densidad de obstáculos intermedia, 1. e., el 
SSP. A nivel de riqueza, el SSP no se diferenció del B, 
esta variable fue mayor en el SC lo que contrasta con 
la literatura que documenta la relación positiva entre 
la complejidad de la estructura vegetal y la riqueza de 
especies (August, 1983). Así, nuestros resultados indi- 
can que la disminución del área boscosa producto de 
la expansión de los SC en la zona de estudio ha tenido 
un efecto negativo sobre las especies propias del bos- 
que. Inventarios de murciélagos del piedemonte de 
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la Cordillera Oriental sugieren que estos ensambles 
son ricos en murciélagos insectívoros (Ramírez-Chaves 
et al., 2013; Sánchez, 2017), lo que es confirmado por 
nuestros resultados. Probablemente esta riqueza, que 
implica redundancias ecológicas (Walker, 1992), per- 
mite que haya especies que hoy, después de la intensa 
perturbación, puedan ocupar los ambientes abiertos 
creados por la ganadería. 


Por otra parte, otro factor que puede influir en los re- 
sultados obtenidos está relacionado con el método de 
muestreo. La naturaleza abierta del SC permite que 
el micrófono capte señales sin interferencias a mayor 
distancia, debido a que hay una menor atenuación de 
los llamados en ambientes con ramas y vegetación 
(Kraker-Castañeda et al., 2013). Adicionalmente, la im- 
posibilidad logística para ubicar el micrófono a más de 
tres metros del suelo dentro del sotobosque probable- 
mente no permitió capturar las señales que emiten los 
murciélagos insectívoros que vuelan por encima del do- 
sel. Sin embargo, en otras áreas del neotrópico, y usan- 
do métodos acústicos, se ha podido establecer que los 
murciélagos insectívoros usan en mayor proporción los 
fragmentos de bosque que las áreas abiertas, por lo que 
estas herramientas parecen ser adecuadas para evaluar 
efectos como los presentados aquí (Estrada et al., 2004; 
Kraker-Castañeda et al., 2013). 


Por otra parte, el uso de detectores ultrasónicos ha per- 
mitido ampliar considerablemente los inventarios de 
murciélagos neotropicales (Azofeifa et al., 2019; Ochoa 
et al., 2000; Zamora-Gutierrez et al., 2016). El listado 
de murciélagos insectívoros aéreos presentado en este 
trabajo es el primero obtenido por métodos acústicos 
para la Orinoquia, y sugiere la presencia de dos nue- 
vos registros para la región: C. greenhalli y P. nasutus 
(Ramirez-Chaves et al., 2016; Solari et al., 2013). Se re- 
quieren exploraciones adicionales, con obtención de lla- 
madas de referencia, y capturas con redes de dosel para 
confirmar estos hallazgos. 


Uso de las coberturas. Los murciélagos insectívoros 
aéreos hicieron un uso diferencial de las coberturas 
evaluadas. El SC fue la cobertura con mayor uso, se- 
guida del SSP y finalmente el bosque. Estos resulta- 
dos concuerdan con otros estudios en áreas altamente 
perturbadas, en las que se ha registrado que la activi- 
dad general es afectada por la composición del paisaje 
(Gerth & Chelsvig, 2008). Sin embargo, nuestros resul- 
tados contrastan con un estudio en México en el que los 
murciélagos insectívoros usaron más los fragmentos de 
bosque que las áreas abiertas (Estrada ef al., 2004). En 
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otro ejemplo de México, los murciélagos insectívoros 
no tuvieron diferencias en el uso de bosques y pasti- 
zales (Kraker-Castañeda et al., 2013). Así, en el futuro 
sería recomendable comparar diferentes sitios con SC, 
SSP y bosques para confirmar si el efecto reconocido en 
este trabajo corresponde a un patrón, o es una excep- 
ción asociada a las particularidades del sitio de estudio. 


Detectamos un efecto asociado a los meses de muestreo 
sobre el total de fases de búsqueda y terminales, pero 
no encontramos una conexión con la estacionalidad de 
la lluvia en la región, excepto que el mes con menor ac- 
tividad fue el periodo seco. Otros estudios también re- 
portan un efecto de la temporalidad en el uso de hábitat 
de los murciélagos insectívoros y se ha argumentado 
que el efecto es explicado por el régimen de lluvias y su 
relación con la disponibilidad de insectos (Barros et al., 
2014; Williams-Guillén & Perfecto, 2011). Por lo tanto, 
estudios de mayor duración con una mejor representa- 
ción de cada fase del ciclo hidrológico, y con mediciones 
de la abundancia de insectos en la zona, pueden ayudar 
a explicar mejor la variación temporal que detectamos. 


Uso de las coberturas por especie. Las seis especies ana- 
lizadas usaron de manera diferencial las tres coberturas. 
A excepción de S. leptura que, en términos de fases ter- 
minales, prefirió usar el SSP sobre el SC y el bosque, el 
uso de las otras cinco especies fue SC > SSP > Bosque. El 
uso de hábitat en los murciélagos insectívoros se relacio- 
na con características particulares en la movilidad y las 
características de los pulsos de ecolocalización (Bader 
et al., 2015; Schnitzler & Kalko, 2001). Así, cada especie 
de murciélago tiene un hábitat en el cual se desempefia 
mejor (Kalko et al., 2008; Marinello & Bernard, 2014), y 
se pueden reconocer murciélagos adaptados a bosques 
y otros murciélagos adaptados a espacios abiertos, y 
cada uno de ellos usa en mayor proporción las cober- 
turas que más le conviene (Estrada-Villegas et al., 2010; 
Estrada et al., 2004; Jung & Kalko, 2010. 


Para la familia Vespertilionidae, las especies con mayor 
actividad en el paisaje ganadero fueron Eptesicus I y M. 
nigricans. A diferencia de las demás especies neotropica- 
les del género Myotis, M. nigricans prefiere los ambientes 
abiertos (Fenton & Bogdanowicz, 2002; Siemer et al., 
2001), lo que concuerda con nuestros resultados. Por 
otra parte, se ha indicado que las especies del género Ep- 
tesicus tienden a tener una carga alar baja, una relación 
de aspecto media y modulaciones en las frecuencias 
emitidas, que les permite tener una alta maniobrabili- 
dad, un vuelo suave y reconocer obstáculos de manera 
eficiente (Bader et al., 2015; Norberg & Rayner, 1987). 
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Estas caracteristicas son apropiadas para espacios de 
borde o cerrados, lo que contrasta con lo encontrado 
en este estudio. Asi, nuestros resultados probablemen- 
te son una evidencia de la versatilidad de los Eptesicus. 
La captura de especimenes de este género en el area de 
estudio debe permitir una descripción más detallada de 
sus características y así ayudar a resolver estas dudas. 


En general, los Emballonuridae neotropicales son ca- 
paces de detectar presas en el borde y claros de los 
bosques (Jung et al., 2014; Schnitzler & Kalko, 2001). 
Sin embargo, P. macrotis tiene una carga alar media y 
una relación de aspecto alta, que le permite volar rá- 
pido y aprovechar espacios abiertos (Estrada-Villegas 
et al., 2010; Marinello & Bernard, 2014; Marques et al., 
2016) lo que concuerda con nuestros resultados. Las 
especies S. leptura y C. brevirostris tienen una carga alar 
baja y una relación de aspecto alta, que les permite te- 
ner una alta maniobrabilidad a vuelos relativamente 
rápidos y aprovechar los bordes de los bosques y los 
claros (Bader et al., 2015; Kalko et al., 2008; Marques 
et al., 2016). A pesar de esto C. brevirostris usó con ma- 
yor frecuencia las zonas abiertas, mientras que S. leptura 
usó más el SSP. Esto podría relacionarse con el hecho 
de que C. brevirostris (machos 7.9 g; hembras 10.15 g) es 
más grande que S. leptura (machos 4.4 g; hembras 5.7 g) 
(Eisenberg, 1989). Adicionalmente, es interesante que 
en el ensamble la única que parece preferir el SSP sea S. 
leptura, y en el futuro debería evaluarse si esto facilita 
su coexistencia con las otras especies. 


El único Molossidae analizado fue M. cf. rufus. Molossus 
rufus tiene una relación de aspecto y carga alar media, 
pulsos de ecolocalización largos y con baja frecuencia, 
y en algunos casos con alternancia, que les permite vo- 
lar rápidamente y ser eficientes en capturar presas en 
zonas semi-abiertas (Bader et al., 2015; Marques et al., 
2016). Esto concuerda con el bajo registro de actividad 
de este murciélago en el bosque. 


De acuerdo a lo discutido, con base en la morfología alar 
y en las señales de ecolocalización, varias de las espe- 
cies de murciélagos insectívoros registrados deberían 
ser capaces de explotar de manera eficiente el SSP, pero 
su actividad se concentra principalmente en el SC. En 
el momento, solo podemos especular, pero las zonas 
de SSP y bosque tienen no solo más obstáculos, sino 
tal vez más sitios de percha para depredadores como 
búhos y halcones que cazan al atardecer o en la noche. 
Así, estudios futuros deberían explorar en conjunto el 
uso del espacio por depredadores de murciélagos que 
podrían explicar el mayor uso de las áreas abiertas 
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En conclusión, la heterogeneidad espacial del paisaje 
ganadero estudiado afectó la estructura del ensamble de 
murciélagos insectívoros aéreos, así como el uso de há- 
bitat de algunas especies. Nuestros resultados resaltan 
la utilidad de los métodos bioacústicos para reconocer 
los efectos de los cambios a nivel de paisaje tanto sobre 
el ensamble, como sobre las respuestas de las especies 
que lo conforman. En la medida en que la industria 
agropecuaria sea más productiva y sostenible, el uso de 
técnicas de monitoreo acústico estandarizado podrá ser 
usado en áreas más extensas para proveer indicadores 
del impacto y de la sostenibilidad de la industria agrí- 
cola (Doohan et al., 2019). 
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Anexo 


Anexo la: El 76.6 % de las identificaciones realizadas para 
la familia Molossidae fueron asignadas correctamente 
de acuerdo a los análisis discriminantes. Dos funciones 
discriminantes explicaron el 97.5 % de la varianza de los 
datos (Anexo 1b). Los parámetros acústicos con mayor 
peso en las dos funciones discriminantes fueron frecuen- 
cia final (Lambda de wilks = 0.006; p < 0.001), frecuencia 
de mayor energía (Lambda de wilks = 0.003; p < 0.001), 
ancho de banda (Lambda de wilks = 0.003; p < 0.001), in- 
tervalo de pulso (Lambda de wilks = 0.003; p < 0.001) y 
frecuencia máxima (Lambda de wilks = 0.002; p < 0.001) 
(Anexo 1c). Para la familia Emballonuridae, las funcio- 
nes discriminantes asignaron correctamente el 93.2 % de 
las identificaciones realizadas. Dos funciones discrimi- 
nantes explicaron el 99. 2% de la varianza de los datos 
(Anexo 1b). Los parámetros acústicos con mayor peso 
en las dos funciones discriminantes fueron: frecuencia 
de mayor energía (Lambda de wilks = 0.002; p = 0.000), 
frecuencia inicial (Lambda de wilks = 0.005; p < 0.001), 
Duración (Lambda de wilks = 0.003; p < 0.001), frecuencia 
máxima (Lambda de wilks = 0.003; p < 0.001), frecuencia 
mínima (Lambda de wilks = 0.003; p < 0.001) e intervalo 


) 
) 
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del pulso (Lambda de wilks = 0.002; p < 0.001) (Anexo 1c). 
Para la familia Vespertilionidae y Thyropteridae, las fun- 
ciones discriminantes asignaron correctamente el 82.3 % 
de las identificaciones realizadas. Dos funciones discri- 
minantes explicaron el 98.1 % de la varianza de los datos 
(Anexo 1b). Los parámetros acústicos con mayor peso en 
las dos funciones discriminantes fueron: Frecuencia ini- 
cial (Lambda de wilks = 0.006; p < 0.001), frecuencia final 
(Lambda de wilks = 0.005; p < 0.001), frecuencia máxima 
(Lambda de wilks = 0.005; p < 0.001), frecuencia mínima 
(Lambda de wilks = 0.006; p < 0.001), Duración (Lambda 
de wilks = 0.007; p < 0.001) e intervalo de pulso (Lambda 
de wilks = 0.005; p < 0.001) (Anexo 1c). Para los sonotipos, 
las funciones discriminantes asignaron correctamente el 
93.3 % de las identificaciones realizadas. Dos funciones 
discriminantes explicaron el 95.4 % de la varianza de los 
datos (Anexo 1b). Los parámetros acústicos con mayor 
peso en las dos funciones discriminantes fueron: frecuen- 
cia máxima (Lambda de wilks = 0.002; p < 0.001), ancho 
de banda (Lambda de wilks = 0.004; p < 0.001), frecuen- 
cia inicial (Lambda de wilks = 0.003; p < 0.001), frecuen- 
cia final (Lambda de wilks = 0.002; p < 0.001) y duración 
(Lambda de wilks = 0.002; p < 0.001) (Anexo 1c). 


Anexo 1b. Parámetros generales de las funciones discrimantes por cada familia de murciélagos registrados por métodos bioa- 


custicos en un paisaje ganadera del piedemonte llanero colombiano. En todos los análisis discrimantes se explicó más del 90 


% de la variación. Var (%) = Porcentaje de variación explicada por función discriminante. Acum (%) = Porcentaje acumulado 


de las funciones discriminantes. 


Annex 1b. General parameters of the discriminating functions for each family of bats recorded using bioacoustic methods in a cattle ranching 


landscape in the Andean foothills of Colombia. In all discriminating analyzes, more than 90% of the variation was explained. Var (%) = Per- 


centage of variation explained by discriminant function. Accum (%) = Accumulated percentage of discriminant functions 


Molossidae 
f.d Autovalor Var (%) Acum. (%) Correlación canónica 

54.8 93.7 93.7 0.991 

2 pu 3.8 DE 0.831 
Emballonuridae 

104.4 98.1 98.1 0.995 

2 1.2 1.1 99.2 0.735 
Vespertilionidae 

50.7 96.0 96.0 0.990 

2 1.1 2:4 98.1 0.726 

Sonotipos 
37.6 81.6 81.6 0.987 
2 64.1 13.9 95.4 0.930 
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Anexo 1c. Parámetro acústicos mejor explicados en las funciones discriminates para cada familia. f.d = función discriminante. 


Annex 1c. Best explained acoustic parameters in the discriminant functions for each family. f.d = discriminant function. 


Molossidae 

Variable f.d 1 d2 

Frecuencia final 0.960 0.121 

Frecuencia Mayor energía 0.807 -0.391 
Frecuencia máxima 0,646 -0.340 
Ancho de banda 0.117 -0.297 
Intervalo de pulso -0.155 0.490 

Emballonuridae 

Frecuencia Mayor energía 0.886 0.212 
Frecuencia final 0.877 -0.464 
Frecuencia máxima 0.834 0.248 
Duración -0.062 0.111 

Frecuencia mínima 0.597 0.289 
Intervalo de pulso -0.095 0.120 

Vespertilionidae 
Frecuencia inicial 0.903 -0.206 
Frecuencia mínima 0.819 -0.086 
Frecuencia máxima 0.531 -0.135 
Intervalo de pulso -0.174 0.574 
Duración -0.250 0.580 
Frecuencia final 0.577 0.129 
Sonotipos 

Frecuencia máxima 0.729 -0.080 
Frecuencia final 0.656 -0.109 
Frecuencia inicial 0.278 -0.512 
Ancho de banda 0.231 0.614 

Duración -0.053 -0.214 
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Anexo 1c. Agrupación de las especies de acuerdo a las funciones discriminantes generadas por familia. funtamos las familias 
Vespertilionidae y Thyropteridae por la similitud en los parámetros acústicos. 


Annex 1c. Species grouping based on the discriminant functions generated for each family. We pooled the data from the families Vesperti- 
lionidae and Thyropteridae due to the similarity in their acoustic parameters. 
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Anexo 2. Numero de fases de búsqueda y de fases terminales de murciélagos insectívoros aéreos registradas en tres tipos 
de coberturas en un paisaje ganadero en el Meta, Colombia. En negrilla las especies elegidas para analizar el uso del espacio 


entre las coberturas evaluadas. Bosque de piedemonte = B, sistema silvopastoril = SSP y sistema convencional = SC. 


Annex 2. Number of search phases and terminal phases of aerial insectivorous bats recorded in three types of covers in a livestock landscape 
in Meta, Colombia. In bold the species chosen to analyze the use of the space between the evaluated covers. Piedmont forest = B, silvopas- 
toral system = SSP and conventional system = SC. 


B SC SSP B SC SSP 


Especie > Total ——_—____——_ Total 
Fases de busqueda Fases terminales 
Cynomops greenhalli 3 3 
Eumops glaucinus 7 52 85 144 3 8 11 
Molossops temminckii 2 2 
Molossus cf. rufus 55 387 521 963 1 136 50 187 
l Molossus cf. molossus 11 86 35 132 4 2 6 
Molossidae Aro eae 
Nyctinomops cf. laticaudatus 
Nyctinomops macrotis 3 2 5 2 2 
Promops nasutus 
cf. Tadarida brasiliensis 11 9 20 
Molossidae A 11 11 1 1 
Eptesicus | (23) 1545. 25615 3127 157 98 132 367 
Lasiurus sp. 4 2 6 
M Myotis nigricans 104 431 187 722 3 19 8 30 
Vespertilionidae ee ee 
Myotis riparius 14 27 41 
Vespertilionidae A 2 2 1 1 
Vespertilionidae B 6 7 13 1 1 
Peropteryx macrotis 51 283 64 398 3 74 8 85 
Cormura brevirostris 21 85 73 179 1 6 li 
ptas Diclidurus ws 9 4 13 1 1 
Saccopteryx bilineata 56 2 2 60 15 3 18 
Saccopteryx leptura 84. 87 298 469 8 13 87 108 
Saccopteryx sp. 2 2 
Thyropteridae Thyroptera sp. 4 4 
Sonotipo 1 1 6 7 
Sonotipo 2 5 5 
Sonotipo 3 2 2 
Sonotipo 4 2 2 
Sonotipo 5 18 18 
Sonotipo 6 2 2 
Sonotipo 7 2 1 3 
No. Id 23 60 22 105 2 2 
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Anexo 3. Porcentaje de la contribución de cada especie de murciélago insectivoro aéreo a la medida de disimilaridad de 
Bray-Curtis entre tres tipos de coberturas en un paisaje ganadero en el Meta, Colombia. Bosque de piedemonte = B, sistema 
silvopastoril = SSP y sistema convencional = SC. Cont. Corresponde al porcentaje de disimilitud contribuido por cada espe- 


cie. Acu. corresponde al porcentaje de disimilitud acumulado. 


Annex 3. Percentage of the contribution of each species of aerial insectivorous bat to the Bray-Curtis dissimilarity measure between three types 
of cover in a livestock landscape in Meta, Colombia. Piedmont forest = B, silvopastoral system = SSP and conventional system =SC. Cont. 
Corresponds to the percentage of dissimilarity contributed by each species. Acu. corresponds to the percentage of dissimilarity accumulated. 


B vs SC B vs SSP SSP vs SP 
Disimilitud = 74.9% Disimilitud = 72.7% Disimilitud = 52.7% 
Especie Cont. Acu. Especie Cont. Acu. Especie Cont. Acu. 
(^) () (^) ^ (Qo) (^) (o) 
Eptesicus I 44.66 | 44.66 X EptesicusI 404 404  EptesicusI 38.7 | 3847 
Molossus cf. rufus 16.0 60.6 Molossus cf. rufus 18.6 59.1 Myotis nigricans 147 53.4 
Myotis nigricans 125 731 Saccopteryx leptura 11. 705 Molossus cf. rufus 15:5 067-2 
Peropteryx macrotis 9.6 827 Myotis nigricans 8.6 79.1 Peropteryx macrotis 9.6 76.8 
Cormura brevirostris 4.3 87.1 Eumops glaucinus 52 84.2 Saccopteryx leptura 8.8 85.6 
Molossus cf. molossus 2.7 89.8 Peropteryx macrotis 4.0 88.3 Cormura brevirostris 52 90.8 
Saccopteryx leptura 2.6 92.3 Saccopteryx bilineata 3.0 91.3 Eumops glaucinus 3.1 93.8 
Eumops glaucinus 2: 94.6 Cormura brevirostris 2.6 93.8 Molossus cf. molossus 2.5 96.3 
Saccopteryx bilineata 2:1 96.7 Molossus cf. molossus — 1.7 95.5 Myotis riparius 1.0 ES 
Sonotipo 5 0.6 97.3 Myotis riparius 1.4 96.9 cf. Tadarida brasiliensis 0.6 97.8 
Myotis riparius 0.6 97.9 Sonotipo 5 1.0 97.9 Vespertilionidae B 0.4 98.3 
cf. Tadarida brasiliensis 0.2 98.1 Vespertilionidae B 0.5 98.3 Diclidurus albus 0.4 98.6 
Molossidae A 0.2 98.4 cf. Tadarida brasiliensis 0.4 98.8  Molossidae A 0.3 98.9 
Diclidurus albus 0.2 98.6 Diclidurus albus 0.3 99.1 Lasiurus sp. 0.2 99.1 
Vespertilionidae B 0.2 98.8 Thyroptera sp. 0.2 99.3  Sonotipo 1 0.2 99 
Sonotipo 1 0.2 98.9 Saccopteryx sp. 0.2 99.5 Nyctinomops macrotis 0.2 99.4 
Thyroptera sp. 0.1 991  Sonotipo 2 0.2 99.6 Saccopteryx bilineata 0.2 99.6 
Sonotipo 2 0.1 99.2 Sonotipo 7 0.1 99.7 | Vespertilionidae A 0.1 99.7 
Lasiurus sp. 0.1 99.3 Lasiurus sp. 0.1 998  Sonotipo 6 0.1 99.7 
Sonotipo 7 0.1 99.4 Nyctinomops macrotis 0.1 99.9 Promops nasutus 0.1 99.8 
Saccopteryx sp. 0.1 99.5  Sonotipo4 0.1 99.9  Sonotipo 7 0.1 99.9 
Nyctinomops macrotis 0.1 99.6 Sonotipo 3 0.1 100.0 Cynomops greenhalli 0.1 99.9 
Vespertilionidae A 0.1 99.7  Sonotipo 1 0.0 100.0 Molossops temminckii — 0.1 100.0 
Sonotipo 6 01 997 Vespertilionidae A 00 100.0 IAN En 0.0 100.0 
Promops nasutus 0.1 99.8 e 2G 0.0 100.0 Thyroptera sp. 0.0 100.0 
Cynomops greenhalli 0.0 99.8  Sonotipo 6 0.0 100.0 Sonotipo 5 0.0 100.0 
Molossops temminckii — 0.0 99.9 — Cynomops greenhalli 0.0 100.0 Sonotipo 4 0.0 100.0 
Sonotipo 4 0.0 99.9 Molossops temminckii 0.0 100.0 Sonotipo 3 0.0 100.0 
Sonotipo 3 0.0 100.0 Promops nasutus 0.0 100.0 Sonotipo 2 0.0 100.0 
ee E 00 . 1000 Molossidae A 0.0 . 1000 Saccopteryx sp. 00 100.0 
aticaudatus 
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Anexo 4. Efecto de las coberturas y los meses de muestreo en el uso de habitat de seis especies de murciélagos insectivoros. 
Todas las especies hicieron un uso diferencial de las coberturas. El factor mes se agregó para controlar el efecto. g.l = grados 
de libertad. d.e = desviación estándar. 


Annex 4. Effect of the covers and the sampling months on the habitat use of six species of insectivorous bats. All species made differential 
use of the covers. The month factor was added to control the effect. g.1 = degrees of freedom. d.e = standard deviation. 


Especie Fase Factores g.l d.e p 
Cobertura 3 24602.040 0.000 
F. búsqueda 
Meses 3 58.680 0.000 
Eptesicus I 
Cobertura 3 1323.478 0.000 
F. terminal wm 
Meses 3 131.697 0.000 
Cobertura 3 3752.740 0.000 
F. busqueda M ————————————— 
Meses 3 268.940 0.000 
M. nigricans 
Cobertura 3 14.985 0.001 
F. terminal ————————————————————————————————————— 
Meses 3 7.333 0.063 
Cobertura 3 1725.471 0.000 
F. búsqueda TTT TFT 
Meses 3 69.496 0.000 
P. macrotis 
Cobertura 3 201.359 0.000 
F. terminal E ———————————————————————————————— 
Meses 3 35.010 0.000 
Cobertura 3 2076.397 0.000 
F. búsqueda —— 
Meses 3 174.241 0.000 
S. leptura 
Cobertura 3 259.868 0.000 
F. terminal > 
Meses 3 103.150 0.000 
Cobertura 3 455.088 0.000 
F. búsqueda ——4—-—-—»————————————————— 
Meses 3 95.351 0.000 
C. brevirostris 
Cobertura 3 26.389 0.000 
F. terminal o 
Meses 3 13.156 0.000 
Cobertura 3 5688.115 0.000 
F. búsqueda E 
Meses 3 144.481 0.000 
M. cf. rufus 
Cobertura 3 623.733 0.000 
F. terminal S 
Meses 3 68.833 0.000 
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